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Abstract: We present a quantitative analysis of the constraints on the total equa- 
tion of state parameter that can be obtained from measuring the red shift evolution 
of the cosmic shear. We compare the constraints that can be obtained from mea- 
surements of the spin two angular multipole moments of the cosmic shear to those 
resulting from the two dimensional and three dimensional power spectra of the cos- 
mic shear. We find that if the multipole moments of the cosmic shear are measured 
accurately enough for a few red shifts the constraints on the dark energy equation of 
state parameter improve significantly compared to those that can be obtained from 
other measurements. 
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1. Introduction 

The observational evidence that the expansion of the universe is being accelerated 
by dark energy is strengthening continuously [l], [2|, |3|, f|]. On the other hand, the 
theoretical efforts to understand the nature of dark energy and the underlying reason 
for the accelerated expansion were less successful so far. It has become clear that 
in order to make further theoretical progress better methods of determining the 
expansion history of the universe are needed. The kinematic distance probes in 
the homogeneous and isotropic universe, such as luminosity distance (g?l), angular 
distance etc., are limited because they rely on light emitted from distant sources, and 
hence measure an integral over the expansion history ||. Consequently even very 
precise future measurements of the kinematic probes will not be sufficient to extract 
precise model-independent information on the expansion history. In addition to the 
kinematic probes one can use dynamical (fluctuations) probes to help determine the 
expansion history of the universe. Cosmological perturbations provide, through their 
dependence on the homogeneous and isotropic background, an independent source 
of information |5|, [7], A promising probe is the weak gravitational lensing induced 
by the perturbations. The three dimensional cosmic shear is a redshift dependent 
tensor that measures the shape distortion of distant galaxies as the light emitted 
from them propagates through the perturbed universe |9|, The two dimensional 



cosmic shear has been measured in weak lensing experiments jH, [E| 13] where the 
best cosmological constraints are obtained from |14|, [15[]. Perliminary studies of 3D 
analysis of current data have been done |16], |7j and planned programs (LSST[|H 



SNAP/JDEM, DUNE/EUCLIDP], GO]) are expected to have some three dimensional 



- 1 - 



capabilities. Several reviews of the theory and observations of weak lensing provide 
an excellent description of the field [^TJ |23| . 

Previously several investigations on the possible use of the three dimensional 
cosmic shear to determine the expansion history of the universe or the nature of dark 



energy were carried out J23], |25], gS[ |27|, |29], |3(], [31], [32| . The standard approach to 
the analysis of properties of dark energy using weak lensing relates the shear to the 
matter power spectrum. Choosing a parametrization for the dark energy equation 
of state one can then make numerical estimations for the cosmological parameters, 
using the results of simulations to take into account the contribution of the non-linear 
regime. Issues concerning the choice of the equation of state parameterizations in 
this context have been discussed in |33, 34]. In order to analyze the three dimensional 
cosmic shear functional dependence on the dark energy equation of state in a model 
independent way one needs observables that can determine in a reliable way the time 
dependence of the perturbations' amplitude and not only their spectrum at a fixed 
redshift. Such observables can have a simple functional dependence on the expansion 
history and hence can provide much more information on its time-dependence than 
the measurement of a single spectrum. If the cosmological perturbations are small, 
the cosmic shear depends linearly on them. In this linear regime both observations 
and theoretical calculations can be done accurately. 

In this paper we analyze the possibility to constrain the total equation of state 
of the universe w to t(z) (the ratio of the total pressure of the universe to its energy 
density) in a model-independent way by cosmic shear measurements. By "model- 
independent" we mean the functional dependence of the cosmic shear on the expan- 
sion history of the universe rather than a specific parametrization. In the context 
of this paper this is equivalent to determining the sensitivity of the angular spec- 
trum of the three dimensional cosmic shear to changes in the total equation of state 
parameter. The question that we wish to pose and answer is: What is the quantita- 
tive precision on model-independent information about the expansion history of the 
universe that can be obtained from cosmic shear measurements with a given accu- 
racy? In particular, we wish to determine whether the angular power spectrum of 
the three dimensional cosmic shear is more sensitive to changes in the total equation 
of state parameter than other statistics, and if it is, then by how much. Our re- 
sults allow us to estimate the accuracy goal needed for shear measurements so they 
can improve on other accurate tests such as luminosity distance measurements or 
CMB measurements. Our approach is mostly relevant when trying to estimate the 
prospects of future lensing surveys for constraining the evolution of the universe in a 
model independent way. Our results suggest that precise measurements of the three 
dimensional cosmic shear at different redshifts will significantly enhance our ability 
to constrain the expansion history of the universe in a model-independent way. 

The paper is organized as follows. In Sect. |2| we review the theoretical material 
from PS[] on which the numerical results, Sect. [3, will be based. In Sect. [| we com- 
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pare the cosmological constraints that can be obtained from the angular multipole 
moments of the three dimensional cosmic shear to those that can be obtained from 
the 2D and 3D spectra of the cosmic shear and Sect. | contains our conclusions. 



2. Theoretical background 

In this section we review the theoretical results from [[JIJ (where many additional 
details can be found) on which the subsequent numerical analysis is based. The three 
dimensional cosmic shear, 7, is related to the metric perturbations, $, through the 
integral expression 

w 

7= dw 55$. 2.1 

J fK{W)f K {w') 


Where w is the radial distance coordinate and fK{w) = K~ 1 / 2 si^K^-^w) , w or 
(_iT)-V2 S inh((-ir) 1 / 2 w) for a closed (K > 0), flat (K = 0) or open (K < 0) 
universe respectably, and 5 is the spin-weight operator . See P5j Section 3 for 
details. 

The metric perturbation is related to the matter density perturbation 5 m by the 
Poisson equation 

V 2 $ = ^^(l + ^ m . (2.2) 

In the linear regime the solutions for 5 m are described by a growth function g(z) 
and the spectrum at some initial given redshift of 5 m (zo). In many analyses, one is 
interested also in the nonlinear regime of 5 m . In this case, one has to use simulations 
or non-linear approximations and interpolate between the linear and non-linear so- 
lutions. We wish to avoid using the non-linear regime as much as possible so we can 
take advantage of the better accuracy of the theoretical calculations in the linear 
regime. Previously, some proposals on how to sidestep the problems of modeling the 
non-linear regime were put forward. It was argued in |27[] that one should elimi- 
nate the dependence on the growth factor and concentrate on geometric quantities. 
A similar approach was followed in pS . In this case, of course, the measurements 



will have similar limitations as the measurements of other kinematic probes (as an 
example one can compare the resulting constraints in [|36| with those in [[321 ). 



Our approach is different. Since it is the perturbations of the metric that shear 
really depends on, we solve directly for $ rather than for 5 m . While the solution 
for the growth function with an arbitrary time dependent equation of state can be 
calculated numerically, the solutions for $ have a simple functional dependence on 
w to t{z). Further, the linear approximation in this case requires that |$| -C 1 but not 
necessarily that 5p m be small. At the epoch of matter-radiation equality the density 
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and metric perturbations are very small, with an amplitude of the order of 10 -5 . 
Through matter domination the metric perturbations are constant while the matter 
density perturbations grow as the scale factor. Today, 5p m can reach and on some 
scales exceed unity. The metric perturbation $, on the other hand, stays frozen at 
a small value. The value of $ is small also for very small scales, for example, in a 
galaxy with 5p m ~ 100 the value of $ is less than 10~ 3 . 

For late times, we have to take into account the effects of dark energy. $ is no 
longer constant, however, the solution for $ is separable |35 



VP 



d7]a 2 (?f)(l + Wtat(rj)). 



(2.3) 



where 77 is the conformal time. In red-shift space the time-dependence of the pertur- 
bation is thus $t(z) which has a simple functional dependence on w to t through the 
first order differential equation 



2(1 + z) 



5 + 3w 



tot\ 



l+t»tot(z) 
5+3wtot(«) ' 



(2.4) 



This solution can be used to calculate the growth function for an arbitrary 
w tot (z). The solution can also be used to construct from eq. fl2.1|) the spin two angular 
power spectrum of the shear 



2 Cl(w 1 ,W 2 ) = (47T 
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where Wi )2 are the distances of the sources and q is the wave number of the pertur- 
bation. For a flat spectrum 



/(?) = A 



2n 



(2.6) 



T q 2 (r]i n ) being the transfer function at early stages of the matter domination epoch 



(denoted by an initial time r] in ) and A is the primordial amplitude. In [^5] we 
have found that the 2 C^s are not sensitive to changes in the shape of the initial 
spectrum. The evolution with red shift of the 2 C\ can be evaluated only if we know 
the distance-redshift relation (written here for a spatially flat universe) 



w 



d L {z) 
l + z 



1+2 



Hn 



dx' exp 



dlnx (1 + w to t(x)) 



(2.7) 



In general, d^z), depends on the total equation of state w to t, the spatial curvature 
K and the value of H . If the luminosity distance is known to a good accuracy (equal 
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to, or better than, that of the shear measurement) then for a fixed angular scale Z, 
2Ci(z) is sensitive to changes only in fl m and w to t(z) making it a good candidate for 
constraining the expansion history of the universe. In what follows we will assume 
that d L (z) is known to about a percent level from supernovae observations. 



3. Estimating the cosmic shear angular multipole moments 

To see if the 2Ci(z) can constrain the expansion history of the universe, several 
questions have to be answered quantitatively. First, we would like to know whether 
by measuring the 2Ci(z) it is possible to improve on the kinematical probes and 
by how much. Then, we would like to compare the 2Ci(z) to other statistics that 
are frequently used in weak lensing analysis. Additionally, since future supernovae 
observations expect an accuracy of 1% in determining the luminosity distance, we 
will see how the constraints on Q m and Wtot from knowing di(z) to one percent 
accuracy compare to those from the red shift evolution of the shear spin two angular 
power spectrum. 

In order to estimate the cosmological constraints we make the following assump- 
tions: That future surveys can measure the 2Ci(z) at accuracies that are cosmic 
variance limited for five different red shift bins in the range 0.2 < z < 1.5 for about 
100 multipoles centered around Z = 100 in the linear regime. Under these conditions 
we expect that the 2C\(z) can be measured at the percent level. 

Since each of the tensor 2 Ci(zys is estimated by 2(1 + 1) independent moments 
as for the case of regular scalar angular moments, the error estimate follows in a 
similar way to the error estimate of any Gaussian field measured over a fraction of 
the sky, f S k y , with a finite number of galaxies (see, for example, |^31 and the appendix 



of Hence, the fractional error for the 2Ci(z) at a given red-shift is given by 



2 Ci(z) y (21 + l)f s ky V n gal (z) 2 Q(z 



2 



A 2 Q(z) 2 f <x 7 



1+ (3-1) 



The important difference here is that eq.( |3.1| ) is a function of red-shift. The shot 



noise term includes the intrinsic variance of the shear of a single galaxy <j* ~ 0.1 
and the source galaxy density per steradian, n ga x(z), as a function of red-shift. We 
assume that the source galaxy density distribution is uniform in angle so the total 
number of galaxies is given by iV ga i = ATrn g3 x(z). 

To estimate the error we will look at the lowest red-shift range where the 2Ci(z) 
are smallest and the source galaxy density is the smallest so the constraints on the 
accuracy will be the strongest. For example, for z ~ 0.2 2^i(z) ~ 0.02 27r/Z 2 
while for z ~ 1.5 iC\(z) ~ 27r/Z 2 |3J|. A standard expected source galaxy density 
has the form z 2 e~^ z ^ Z0 ^ 2 whose median red-shift is somewhat larger than zq. To get 
the shot noise down to a comparable level to the cosmic variance error for z = 0.2 
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one therefore needs iV ga i(0.2) ~ 10 I 2 for each angular multipole. For example, the 
cosmic variance error of a survey covering half of the sky at Z = 100 is ~ 14% and 
one needs about 10 5 galaxies to get the shot noise down to this level. To improve 
on the accuracy of cosmic variance limited measurement one can use binning of 
multipoles so that the number of independent measurements ~ (Imax ~ ^min)- To 
reach an accuracy on the order of 1%, about 100 multipoles around I ~ 100 must 
therefore be measured. To reduce the shot noise to the percent level one needs in 
this case N(0.2) ~ 10 (4o*/ 3 - l miJ 3 ) ~ 107 galaxies. If the distribution of source 
galaxies is as given above with median redshift of about 1.25, one would need about 
2 x 10 8 galaxies. This would be equivalent to a total galaxy source density of about 
2 galaxies per squared arcminute. 

Whether or not future cosmological observations of the three dimensional cosmic 
shear and of luminosity distance will achieve the percent accuracy goals that we have 
described will be determined by future technological advancements in the different 
fields of observational cosmology. However, at the moment, there does not seem 
to be any fundamental theoretical or technological reason that indicates that such 
accuracy goals can not be achieved. 

In view of the above, we let the accuracy with which the 2Ci(z) are assumed 
to be measured to vary between 1% and 10% in our calculations. In addition, as 
previously mentioned, we assume a prior of 1% on the redshift-distance relation (i.e. 
luminosity distance). Since we wish to expose the difference between the sensitivities 
of 2Ci( z ) and cIl(z), we put strong priors (1% error) on K and H when calculating 
the sensitivity of cLl{z) which enhances the sensitivity of di{z). If we select weak 
priors the advantage of 2Ci(z) becomes more pronaunced. 

The expression for the 2Ci{ functional of the total equation of state 

parameter given in Eqs. (|2.4j ), (|2.5|) , does not assume a particular form of the time 
dependence of w to t or wde- However, for the purpose of the numerical comparison 
between the different measures of the cosmic shear we choose, for the moment, a 
specific parametrization. We use Linder's parameterization wde = w + Wi(z/1 + 
z), because it is commonly used in the literature and hence its use facilitates the 
comparison of our results to the results of other theoretical investigations and to 
most of the existing observational results. That said, we wish to emphasize that our 
results can be as easily adapted to any other parametrization. In fact, quite a few 
papers suggest other (improved) parameterizations of wde that are just as easy to 
analyze for the measurements of 2Ci(z). We use the standard likelihood analysis to 
estimate the constraints on w , Wi and Q m . For our fiducial model we choose ACDM 
with Q m = 0.3 The results can be seen in Figure (|l|). As we can see, for the 2C\{z) 
to be competitive the measurement accuracy should be 3% or better. For errors of 
10% the constraints from the shear are weaker than those of the luminosity distance, 
on the other hand when the error is 1% the constraints improve significantly. 

A second important aspect of the behavior of the 2C\{z) is the improvement of 
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Figure 1: The la contours for 2Ci( z ) at I = 100 with assumed errors of 1%, 3% and 10% 
(purple lines). The wider filled contours (green) show the la and 2a likelihood contours 
for supernovae di(z) measurements with 1% accuracy. 

the constraints for higher multipoles. In Fig. [|we can see the error in the equation of 
state parameters decreases as a function of I. The reason for that is the fact that the 
integrand of ( |2.5| ) samples $ at lower redshifts for higher orders of the Bessel function 
and since $ is more sensitive to changes in wtot at the lower redshifts there is more 
sensitivity at higher Vs. This is particularly significant for the w\ parameter which 
determines the time variation of the equation of state in Linder's parametrization. 
This behavior is independent of the fact that for higher multipoles the error is smaller 
and demands only that we use the linear regime. 

The range of multipoles that are linear is larger than can be used in the case of 
the standard two dimensional weak lensing observations. The reason is as follows. 
The matter power spectrum is accepted as linear up to k = O.lMpc -1 which for the 
two dimensional shear power spectrum is about / = 200, above which the standard 
practice is to use non- linear models to estimate it. The metric perturbation $, on 
the other hand, is linear at smaller scales (larger Z's). For wavelengths well inside the 
horizon 5p m / p m is larger than $ by a factor of the order of the square of the ratio 
of the size of the horizon to the wavelength, op m /p m ~ (q/H) 2 &. The larger linear 
range allows us to use the higher multipoles effectively. This is an advantage since 
as can be seen from Fig. (|2]), the higher multipoles are more sensitive to changes 
in the equation of state and going beyond / = 200 improves our constraints. The 
transition to non-linear regime for the matter density is on scales of about 8 Mpc, 




Figure 2: The magnitude of the error (la) in wq (bottom line) and w\ (top) as a function 
oil. 

slightly larger than the scale of clusters. The value for $ on the other hand is much 
smaller even down to scales of 200Kpc, the size of a galaxy. When calculating the 
iCi{z) the smallest scale that considered here is for the case / = 500 (the largest I 
considered in the paper) at red shift of z — 0.1 (the smallest redshift considered in 
the paper). In this case the relevant scale is about 1.2 Mpc which is still extremely 
linear in $. As an example, the expected corrections for $ from non-linear terms at 
k = 0.8Mpc _1 are of the order of 10 -7 . 

4. Shear power spectrum comparisons 

There are many ways to quantify and measure shear, each having its own merits and 
drawbacks. Several two dimensional measures of shear have been measured to date 
and a substantial amount of work on the three dimensional power spectrum has been 
published. The cosmic shear depends on the metric perturbations and thus like the 
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matter density perturbations is a good source for information on <t 8 , Vt m , w to t and 
other cosmological parameters. As previously explained, in the context of cosmologi- 
cal models and explaining the expansion history of the universe, the important point 
for us is how well it is possible to constrain the total equation of state as a function of 
red-shift. From this narrow point of view we ask whether different measures of shear 
perform differently. We will show that the 2^i(z) is the most sensitive to w tat (z) 
and gives the best constraints on it. When planning an experiment one questions 
the accuracy with which any quantity can be measured but we would like to ask a 
different question: Given an accuracy, what constraints can we obtain on the rele- 
vant parameters? There are many different considerations and technical issues that 
take part in evaluating future measurements and the expected constraints they will 
give. The answers vary between different surveys and therefore we take this different 
approach where we try and keep the comparison similar so that the results reflect 
the sensitivity of the measures and not the quality of the experiments. To that end 
we choose an accuracy goal which is (almost) optimal so that the results will reflect 
the best future prospects. 

In this section we present a numerical analysis of the sensitivity of different shear 
measures to changes in the total equation of state. We set an a priori accuracy level 
of 1% (which is the goal accuracy for future cosmological observations) and analyze 
the likelihood for Q m , wo and w± (as in the previous section). The measures we 
choose are common to all the different methods and differ between them in their 
intrinsic functional dependence on w to t{z). The first is the two dimensional angular 
power spectrum [ f22| 

D m mfal f w (w) f i \ 



4c 4 J a?{w) 



W(w) = / dw'G{w') 



f K (w' - w) 
f K (W) ' 

w 

where G(w) is the normalized source distance distribution function. The second is 



the three dimensional power spectrum 10 



CP{k u h) = ^r|^| J k'dkhih, fc)/,^, k) (4.2) 

o 

oo r 

hik^^h J drrj^hr) J dr' (^—L^ 3l {kr')^P^{k-r'), 
o o 

and the third is the zCiiz). 

The assumptions concerning the distance red-shift relation and all the priors are 
similar to those in the previous section. In addition, these quantities depend on the 
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Figure 3: The la contours of 2Ci(z) (purple), C^ 77 (fei, k%) (green) and P K {1) (orange). 

normalization of the power spectrum. Although it is common to use erg, the metric 
perturbation is proportional to the primordial amplitude of perturbations, A. The 
value of A is obtained from CMB measurements to 5% > accuracy. The 2Ci(z) are 
insenstive to the value of A and assuming a prior of 5% then marginilizing over it 
makes no differnce to the wo,wi likelihood. For the 2D and 3D spectrums on the 
other hand margenilizing over A increases the contours. Therfore, in the following 
calculations we keep the value of A fixed, so that it is the difference in sensitivity to 
w to t that is demonstrated. 

We can see from the contours in FigfJ that the iCi^z) is more sensitive to 
changes in w to t and Q m , the reason for this is clear upon examination of the func- 
tional dependence of the different measures on these parameters. The dependence 
of perturbations on the total equation of state comes through the growth function. 
The two dimensional power spectrum is a projection, therefore it contains a double 
integral over §t(z) and the result is degenerate with respect to w tot . This is solved 
by measuring distance (red shift) dependent quantities such as the three dimen- 
sional spectrum. But, although the full spectrum measures best all the cosmo logical 
parameters of the model it measures the shape and evolution of the spectrum simul- 
taneously and therefore is less effective in distinguishing between the information 
from the growth function and from the spectrum's shape. The 2Ci(z) on the other 
hand is insensitive to the k dependence of the spectrum and therefore can not deter- 
mine as directly, but it has a straightforward functional dependence on w to t and fl m 
and therefore gives better constraints on them. 
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For the three dimensional methods we have assumed that a single (binned) multi- 
pole is measured but both quantities can be measured at several (binned) multipoles 
and with different correlations (C^ 1 (k 1 ,k 2 ) , fci ^ k 2 or 2 Ci(zi, z 2 ) , z\ ^ z 2 ) thus 
increasing the final accuracy. This is important since for a given survey magnitude 
the error on P K (l) should be smaller than on 2 Ci(z) because of the higher num- 
ber density of galaxies, but as the number of redshifts grows the number of cross 
correlations grows too and with it the accuracy. 

Multi parameter modeling of w^e 

The main advantage of using the 2 Ci(z) is their explicit and fully functional depen- 
dence on w to t without prior assumptions on the functional form of w tot . The major 
improvement of using them compared to other approaches is in their ability to place 
constraints on the evolution of wtot{z)- Since fundamental theory offers little insight 
as to the functional form of the dark energy equation of state, and to demonstrate 
the sensitivity of the 2 C[(z) to changes in it, we choose a multi dimensional param- 
eterization to further explore the difference between the various measures that we 



have mentioned above. We choose a parametrization of wde IhH EH i n which the 
dark energy equation of state is parametrized by a piecewise constant function. The 
constants Wi are the constant values of w(a) in bins of Aa. We choose (for practical 
reasons) to divide the range 1 < a < 0.2 into 5 bins of width Aa = 0.16, so that in 
redshift space Z{ + 1 = 1/eij we have 

4 

w D e{z) = ^2wiT(zi,z i+1 ). (4.3) 

i=0 

The top-hat function T(zi,z i+1 ) is equal to unity for z { > z > and vanishes 
otherwise. In Fig.(§) we present the la errors for each of the parameters marginal- 
ized over Q m . The calculations where done under the same assumptions as in the 
previous sections. This is a somewhat simpler analysis than that of the "principal 
components" but it is sufficient to demonstrate the relative difference between 



the three measures of cosmic shear. We can see once again that in the first 4 bins the 
2 Ci(z) is significantly more sensitive to changes in wde{z) than the other methods 
mentioned above. We can also see once again that the luminosity distance, even when 
determined to one percent, cannot sufficiently constrain the equation of state. In the 
fifth bin, gives the value for z > 1.7 which is beyond the assumed measured range 
(5 measurements in the range 0.2 < z < 1.5) and therefore the constraints are very 



weak. Comparing our results to existing multi-parameter analyses [37], |39|, [MJ, [41], [42 
the most pronounced difference is the low red-shift sensitivity. A possible origin of 
this difference is that the expression of the angular moments contain the square of 



the galaxy distribution function (see, for example, |29] eqs. (2), (3)) which vanishes 



for low redshifts. This issue should be better understood and we leave it for a future 
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Figure 4: The la magnitude of the five redshift bin values of wde( z ) marginalized over 
Q m for 2Ci(z) (purple), CJ 1 (k±, k%) (green) and P R (l) (orange). The empty bars represent 
the la magnitude for di(z) assuming it is measured at 1% accuracy. The redshift bins are 
delimited by {0, .19, .47, .92, 1.78, 4.0}. 



investigation. For higher redshifts, the value of w is constrained by the 2Ci(z) bet- 
ter than the other statistics and the relative advantage of the 2Ci(z) increases with 
redshift. One may wonder why is it that the 2Ci(z) do better than C^ 1 {k\, k 2 ) since 
they contain essentially the same information. The reason is that the redshift evo- 
lution of the 2Ci(z) measures the growth function directly while which C^(ki, h-2) 
does not. Measuring the growth function is what gives us the better constraints on 
the equation of state of the dark energy specifically. Similar effects can be obtained 
by cross correlating P K (l) at different redshift bins. 

We can observe that the 2 C«(< 2 )'s are indeed more sensitive to a deviation of wde 
from —1 as in the case of a cosmological constant and that the extra sensitivity is more 
significant at the higher redshifts. Thus, the three dimensional shear measurements, 
and specifically the red shift evolution of the spin two angular power spectrum, are 
an important tool in determining the evolution of the equation of state. 
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5. Conclusions 



In this paper (as well as its predecessor p5] ) we have clarified the dependence of 
the cosmic shear on the expansion rate of the universe and have made a preliminary 
quantitative analysis of its sensitivity to changes in the dark energy equation of 
state. Theoretical arguments suggest that the red shift evolution of shear spin weight 
two angular power spectrum multipole moments 2^(2) 's, are the best statistics 
for this purpose. They have the simplest dependence on the metric perturbations 
$, and provide a more direct connection between the measured quantities and the 
theoretically interesting quantities. More specifically, theoretical arguments suggest 
that precise knowledge of the 2 Cj(z)' s should allow to obtain better constraints on 
w to t{z) and the expansion history of the universe. We have demonstrated the validity 
of this suggestion by a preliminary numerical analysis. To evaluate the magnitude 
of improvement that we should expect from 2Ci(z) measurements on future surveys 
a full scale numerical analysis of currently planned programs should be carried out. 

Comparing the constraints that can be obtained by using the 2Ci(z) and other 
measures of cosmic shear, such as the 2D and 3D power spectrum as well as kinematic 
distance probes such as d^z), we see that even when all other parameters are kept 
fixed the 2^(2:) 's are more sensitive to changes in and therefore better for 
putting constraints on Wt t{z). Furthermore, we have shown that up to the maximal 
red-shift at which 2C\{z) is measured, it is possible to get a significant improvement 
on the constraints of iftot- This fact is important when one is trying to determine 
whether the dark energy is a cosmological constant or some other dynamical source. 
The constraints on the current value of wde{Q) = —1 are fairly tight so the expected 
improvement from future observations in the Figure of Merit, for example, will come 
mostly from better constraints on the time evolution of the dark energy equation of 
state parameter. Therefore, the best prospect for deciding whether the dark energy 
is a cosmological constant or dynamical lies in better constraints on wtot{z) at higher 
redshifts, such as those that can be obtained from measuring the 2Cz( 2 0' s - 
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